Three dimensional simulation of fluid flow in X‐ray CT images of porous media     

A. Al‐Omari  E. Masad 《国际地质力学数值与分析法杂志》2004,28(13):1327-1360

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction     

M. Zeghal  U. El Shamy 《国际地质力学数值与分析法杂志》2004,28(14):1361-1383

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

广义相对论中静态荷电球体的内解     

王行翔 《中国天文和天体物理学报》1988,(3)

本文在假设荷电球体内部物质密度为ρ_m=μγ~α,电荷密度为ρ_e=ρ_0γ~(b_e-λ/2)的情况下,严格求解Einstein-Maxwell场方程,求得了静态荷电球体的一个较为普遍的内部解。这个解可以看作是Wyman得到的内部解推广到荷电情况下的结果,并且将Wilson,李鉴增,Santos,Pant和Sah等人给出的许多结果作为自己的特例包括在内。  相似文献   

Stochastic Hill's equations for the study of errant rocket burns in orbit     

Martin Ostoja-Starzewski  James M. Longuski 《Celestial Mechanics and Dynamical Astronomy》1992,54(4):295-303

The problem of errant rocket burns in low Earth orbit is of growing interest, especially in the area of safety analysis of nuclear powered spacecraft. The development of stochastic Hill's equations provides a rigorous mathematical tool for the study of such errant rocket maneuvers. These equations are analyzed within the context of a theory of linear dynamical systems driven by a random white noise. It is established that the trajectories of an errant rocket are realizations of a Gauss-Markov process, whose mean vector is given by the solution of a deterministic rocket problem. The time-dependent covariance matrix of the process is derived in an explicit form.  相似文献   

双三次数值模式与理想Lamb波、龙卷扰动对比模拟分析     

辜旭赞  张兵 《气象科技》2007,35(5):613-620

通过"双三次曲面拟合—时间步长积分—双三次曲面拟合……",实现成为一种双三次曲面拟合新算法数值模式。它对模式大气作非线性(空间二阶可导)描述,是用数值分析(三次样条/双三次曲面)与计算方法近似求解大气运动非线性偏微分方程组。用Z坐标系高分辨率双三次数值模式,采用原始大气运动方程组和采用欧拉时间积分方案,设置理想扰动,对大气声波Lamb波做数值模拟;采用大气不可压假定,除去声波,对大气重力波龙卷扰动做数值模拟。对比模拟发现:理想Lamb波以声速在水平方向传播,扰动能量频散快,仅能在原地持续数秒钟;理想龙卷扰动则在原地可持续约1 min,且波动内圈(外圈)自下至上始终处于风场辐合(辐散)状态,具冷心结构和气压驻波特征,可因水汽进入扰动与借助凝结潜热释放而得到发展。  相似文献   

Estimation of parameters in petrologic materials balance equations     

N. H. Gray 《Mathematical Geology》1973,5(3):225-236

A rigorous treatment of the least-squares estimation of the parameters in petrologic materials balance equations is developed to take into account the uncertainties inherent in the chemical analyses. The choice of the optimal estimation procedure for a particular problem is dependent upon the extent of the petrogenetic understanding of the exact nature of the materials balance involved. Some of the complications entailed in obtaining least-squares estimates are illustrated by the examples of magmatic differentiation and metamorphic reactions.  相似文献   

    

Michael W. D. Davis  Michel David  Jean-Marc Belisle 《Mathematical Geology》1978,10(4):369-374

In mineral deposit evaluation, the estimation of block values can involve the solution of as many as 100,000 systems of linear equations. To minimize computing costs it is important that the program to solve the linear systems be as efficient as possible. Standard library routines tend to contain very general stable methods using pivoting and iteration to ensure high accuracy solutions. Here, a particular method has been developed that gives a substantial improvement in run time at the occasional expense of accuracy and stability. The method is not new—just a simple application of Gaussian elimination in the symmetric case.  相似文献   

三角形网格下求解二维浅水方程的和谐Godunov格式   总被引：13，自引：3，他引：13       下载免费PDF全文

潘存鸿 《水科学进展》2007,18(2):204-209

为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。  相似文献   

大欧拉角的空间直角坐标转换方法探讨   总被引：3，自引：0，他引：3  

胡亚江  杨晓梅  沙月进 《现代测绘》2006,29(6):10-12

首先指出建立使用大旋转角的空间直角坐标转换模型的必要性，接着推导出以罗德里格矩阵为旋转矩阵的大欧拉角空间直角坐标转换的数学模型，然后以布尔沙模型为例、以Visual C#．NET2003编制了计算软件。该算法具有性能稳定、计算收敛速度快、精确度高、转换参数初始值可以为任意值的优点。  相似文献   

地震波的场方程矩阵和能量的正定二次型及其意义   总被引：2，自引：0，他引：2       下载免费PDF全文

牛滨华  孙晟  孙春岩  黄新武  李佳 《地球物理学进展》2007,22(2):353-358

场方程是能够表述半空间地震波场的整体特征.波动方程、速度方程和能量方程.通过分析可知每个场方程都具有各自的“场方程矩阵”.能量方程能够对所有场方程矩阵进行综合和贯通，给出了能量方程以“弹性矩阵”为核心的普适性表达形式.最后，运用矩阵的正定二次型理论阐述了“能量矩阵与弹性矩阵”之间一致的对称性和正定性.能量矩阵蕴含的动态力的平衡关系、速度的时间_空间分布和能量的传播及变化的物理意义，能够从能量矩阵的正定二次型特性表述出来.本文研究分析问题的方法完全适用于复杂介质模型，相关的认识和结论可以拓展到均匀黏弹性各向同性介质、均匀弹性各向异性介质、均匀黏弹性各向异性介质以及比奥饱和流体介质.  相似文献